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Abstract 

On a disk deformed to a non-axisymmetric form, a set of oscillations can be excited 
by their resonant interaction through the disk deformation (Kato et al. 2011). This 
resonant instability process has been proposed to suggest a possible cause of the 
high-frequency quasi-periodic oscillations (HF QPOs) observed in black-hole low- 
mass X-ray binaries. In the present paper, we examine whether the above-mentioned 
wave- wave resonant process can describe the tidal instability and superhump in dwarf 
novae. The results show that the process seems to well describe the observations. If 
this process is really the cause of the tidal instability and superhump, a two-armed 
oscillation with high frequency roughly on the magnitude of three times the orbital 
frequency is present on disks, although its expected amplitude may be small. 

Key words: accretion, accrection disks — dwarf novae — oscillations — reso- 
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1. Introduction 

The origins of superoutbursts and superhumps observed in dwarf novae are now well 
understood by the so-called tidal instability and precession of the eccentric disk deformation 
(Osaki 1985, Whitehurst 1988a,b; Hirose & Osaki 1990; Lubow 1991a,b, 1992, 1994: for review 
see Osaki 1996). The eccentric disk deformation and its precession were found by numerical 
simulations by Whitehurst (1988a,b). The origin of the disk deformation was pointed out to be 
due to a parametric resonance between disk rotation and binary revolution from a test-particle 
approximation (Hirose & Osaki 1990) and later in the framework of hydrodynamics (Lubow 
1991a). The precession of the deformation is found to be due to an one-armed global oscillation 
(Osaki 1985; Hirose & Osaki 1993). 

In a different field of astrophysics, i.e., in the field of black-hole low-mass X-ray binaries 
(BH LMXBs), we know that high-frequency quasi-periodic oscillations (HF QPOs) whose fre- 
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quencies are in the range of 100 to 450 Hz are observed in some sources (e.g., for review, see 
van der Klis 2004 and Remillard & McClintock 2006). The origin of these HF QPOs has been 
studied extensively, since clarification of the origin will give a powerful tool to know the inner- 
most structure of relativistic accretion disks as well as the spin of the central black hole sources. 
Inspite of many efforts there is still no consensus on the origin of HF QPOs. However, one 
of possible models of HF QPOs is excitation of a set of oscillations by their resonant coupling 
through disk deformation (Kato et al. 2011). In this model a deformation of unperturbed disks 
from an axisymmetric state is essential, and it is a kind of catalizer for excitation of oscillations. 

The above-mentioned process of excitation of disk oscillations is rather general. Hence, 
a natural question is whether the tidal instability and superhump in dwarf novae can be in- 
terpreted as a result of excitation of disk oscillations on tidally deformed disks by the above- 
mentioned wave-wave resonant process. The purpose of this paper is to demonstrate this 
possibility and to obtain some hints on refining the models of the HF QPOs. 

2. Outline of Wave- Wave Resonant Excitation Process through Disk Deformation 

Let us assume that oscillations in disks can be decomposed into normal modes. The time 
and angular dependences of the displacement vector, $(r,t), associated with the oscillations 
are factorized as $,(r,t) = £(r, z)exp[i(ut — mip)]. Here, r is the cylindrical coordinates (r,ip,z), 
whose center is at the disk center and the z-axis is the axis of the rotating axis of the disk. 
We consider two oscillations. The set of frequency and azimuthal wavenumber, i.e., (u, m), 
of each oscillation is denoted (lux, mi) and (u 2 , m 2 ). Furthermore, we assume that the disk is 
deformed to non-axisymmetric state with azimuthal wavenumber and the pattern rotates 
with frequency wd, i.e., the set of (cu, m) of the disk deformation is (wd, ^d)- In the case of 
tidal deformation, wd and m D are related by 



where Q* rb is the orbital frequency of the secondary star around the primary star, observed from 
the primary star. To avoid unnecessary complication, m-Q is taken to be a possitive integer, 
i.e., m D = 1,2,3.... 

The above two oscillations can non-linearly interact through the disk deformation, if the 
following resonant conditions are satisfied, i.e., 



where the sign of + or — is possible. Kato et al. (2011) showed that the above two oscillations 
grow simultaneously if both oscillations satisfying conditions fl2]) overlap in their propagation 
regions, and have opposite signs of E/u, i.e., (E\/uj\)(E2/uj2) < 0. Here, E is the wave energy 
defined by, e.g., 



= m D Q* 



orb' 



(1) 



(jJ 2 = Wi±u;d an d m 2 = m 1 ±m D , 



(2) 
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Fig. 1. Schcmatical diagram showing resonant amplification of two oscillations with opposite signs of 
lu — mf2, i.e., (cji — mifl)(uj2 — ?ri2^) < through disk deformation. The additional necessary conditions 
for resonance amplification arc U2 = u)\ ± ud, mi = mi ± mo; and ni — n\ ± tid- 



1,2, 



(3) 



where Uo{r) is the velocity on the unperturbed disk, i.e., uq(v) = (0,rf2(r),0) and the asterisk 
shows the complex conjugate. The condition, (E\/uj\)(E2/uj2) < 0, is roughly equal to (uj\ — 
<O0 This resonant amplification process is shown schematically in figure 1. 

This wave amplification comes from an important general characteristic of coupling 
terms between two oscillations through disk deformation. The characteristic is the commutative 
relations given by equation (3) of Kato et al. (2011) or equation (83) of Kato (2008). Explicit 
expressions for the coupling terms are given by equations (5), (6) and (82) of Kato (2008). 

It is noted here that for the resonance to be realized an additional condition is necessary. 
It is a relation among node mumbers, say n's, of oscillations in the vertical direction. This is 
related to the fact that in the normal mode oscillations, oscillations with different node numbers 
are orthogonal with a certain weighting function .§ In the present paper, however, we are only 
interested in the fundamental mode of oscillations in the vertical direction, i.e., ri\ — ri2 = 0. 
Furthermore, the tidal wave is assumed to have no node in the vertical direction, i.e., n-^ = 0. 
Since all n's are zero, such an additional condition is automatically satisfied and no additional 
condition is necessary here. 



In Kato et al. (2011), instead of the condition, (E\ / 'wi)(i?2 / '^2) < 0, opposite signs of wave energies of two 
oscillations, i.e., E\Ei < 0, is sometimes emphasized as the resonant instability condition, since they are 
interested in oscillations with positive frequencies, i.e., u>i > and u>2 > 0. In some cases of disk deformations, 
like the case of the present paper, one of wi and W2 becomes negative. Hence, (uji — m\^l){u2 — m2^l) < is 
more general than E\E-x < 0. 

For vertically isothermal disks see, for example, Okazaki et al. (1987) and for vertically polytropic disks 
see Silbergleit et al. (2001). 
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3. Application to Tidally Deformed Disks 



3.1. Tidally Deformed Disks 

In the case of dwarf novae, the disk of the primary (central) star is deformed by tidal 
force of a secondary star. The time-averaged part of the tidal potential, $, is given by 

*« = -lr 2 . <*> 

where r is the radial distance from the central star, D is the binary separation, and M s is the 
mass of the secondary star. If the pressure force is neglected, the angular velocity of rotation 
of the disk gas, fi(r), is given by 

fir = — 

where M is the mass of the central star and q = M s /M. Since we are interested in the case 
where q is smaller than unity and in the radial region of r/D < 1 , we have approximately 

n = nK(i-70-£)> ( 6 ) 



/ 1 r 3 \ 



CD*, 

where fix is the Keplerian angular velocity of rotation given by fi K (r) = (GM/r 3 ) 1 / 2 . In this 
disk, the epicyclic frequency defined by 

K 2 = 2fi(^2fi + r— J (7) 

is approximately given by 

K = n K (l-q^. (8) 

In addition to the above-mentioned axisymmetric deformation from the Keplerian one, 
the disk of the central star is deformed by non- axisymmetric tidal waves. The set of (cud, ^d) 
of the deformation is given by (m D fi* rb , m D ). 

3.2. One- Armed Low- Frequency Oscillation with u> — mQ < 

Since the disk rotation is slightly deviated from the Keplerian one by the tidal force [see 
eq. OH])], an one-armed (m = 1) p-mode oscillation on the disk is global and has a low frequency 
(Osaki 1985, see also Kato 1983). The dispersion relation for local oscillations [see, e.g., Kato 
2001 and Kato et al. 2008] shows that the radial propagation region of such low frequency 
oscillation [the frequency and azimuthal wavenumber of the oscillation are are denoted u\ and 
mi(= 1), respectively] is specified by 

uji<VL-k. (9) 

Since fl — k is given by [see equations (0) and OH])] 

Q-K=\qnJ^X, (10) 



A \DJ 



and increases outwards, the propagation region of the oscillation with frequency u)\ is bounded 
inside. The inner edge, r c , of the propagation region where u\ = Q — k is given by 



The outer edge of the propagation region is the outer edge of the disk, r t , which will be specified 
as the radius where the disk is truncated by wave-wave resonant instability. 

The low-frequency one-armed oscillation with frequency u\ is thus trapped in the radial 
region of r c < r < r t . The region is schematically shown in figure 2. It is noted that in this 
propagation region, oj\ — mifi is negative, i.e., lo\ — m\Q. < 0. 

3. 3. High- Frequency Oscillation with uo — mQ > and Efficiency of Coupling 

As the counterpart of the Wi-oscillation described above, we consider here an oscillation 
with uj — mQt > which satisfies the resonance conditions (J2J). The oscillation is taken to be 
a p-mode. Its frequency and azimuthal wavenumber are denoted, respectively, w 2 and m-z, 
and are determined later. A p-mode oscillation with given U2 and m 2 has two propagation 
regions of w 2 > m 2 f2 + k and 002 < m 2 f2 — k (e.g., Kato 2001, Kato et al. 2008). In the former 
propagation region we have a; 2 — m 2 f2 > 0, while in the latter we have U2 — m 2 f2 < 0. Hence a 
p-mode oscillation in the former region is our concern here. 

As will be found later, the resonant instability occurs for oscillations with m 2 = —2. 
Thus, 777-2^2 + k < 0, and as is shown schmatically in figure 2, the cu2-oscillation has U2 < and 
is trapped between the inner edge, r s , of the disk and the radius (the Lindblad resonance) 
specified by 



where the subscript L denotes the value at r = r^. Outside of the oscillations are spatially 
damped. For trapping to occur, a relation among r s , rt,, and o; 2 is necessary as a trapping 
condition. In the present problem, however, we need not to pay particular attention on such 
a condition by the following reason. Near to r s , the Co>2-oscillation has very short wavelength 
in the radial direction, since the difference between cu 2 and (m 2 f2 + k)l becomes large there 
and the difference must be compensated by the radial wavelength becoming short (consider 
the dispersion relation of p-mode oscillations.) In other words, the oscillation with o; 2 is one of 
high ovetones and near to frequency U2 there are many eigen-frequencies densely. This means 
that in practice any U2 can become a eigen-frequency and no relation between r s , rt, and o; 2 is 
unnecessary here. 

For the resonant coupling between two oscillations with (cut, mi) and (U2, 1712) through 
disk deformation to occur, the resonant conditions (j5J) are necessary as mentioned before. 
In addition, for the coupling to be efficient, the position of rt must be in the propagation 
region of the Wi-oscillation, since the coupling efficiency is determined by volume integrations 
of some linear products of displacement vectors £ 1; £ 2 , and £ D , they being displacement vectors 




(11) 



u 2 = m 2 Vtt + «l, 



(12) 
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associated with the 00 1- and ^-oscillations and with the uj-q disk deformation (Kato et al. 
2011). For the coupling terms to become large, i) the region where £ 1 has a large amplitude 
and that where £ 2 has a large amplitude must overlap, and ii) both £ x and £ 2 in the overlapped 
region do not vary in the radial direction with short wavelength. The latter requirement comes 
from the fact that if £ 2 , f° r example, changes in the radial direction with short wavelength a 
volume integration of products among £ 1; £ 2 , an d £d becomes small by cancellation. The above 
consideration concerning the coupling efficiency suggests that the resonant couplng occurs most 
strongly in the case of tl ~ r t by the following reasons. The amplitude of the low-frequency 
one-armd oscillation is large around r = r t (see section 5 and also figure 3 by Hirose & Osaki 
1993). Hence, if > r t , the o; 2 -oscillation has short wavelength around r ~ r t and the coupling 
term resulting from the volume integration becomes small. On the other hand, if tl < r t , 
the amplitude of the o; 2 -oscillation is small around r ~ r t , since the region around r t is the 
evanescent region of the o; 2 -oscillation. Based on these considerations, we adopt = r t and 
takfl 



cu 2 = m 2 fi t + /t t , (13) 
where the subscript t denotes the value at r = r t . 

4. Equation Describing Tidal Truncation Radius r t 

As resonant conditions which lead to realistic cases, we adopt 

Co> 2 = cdi — Wd and m 2 = mi — mrj . (14) 

Then, substituting u>i = (f2 — k) c , m 1 = 1, wd = m D f2* rb , and equation (fTBl) into the first relation 
of equations ( fl4|) . we have as the condition of wave- wave resonant instability 

(m D - 2)n t = m D n* orh -{n- «) t - (n - «) c , (15) 

where the subscript c denotes the value at r = r c , and (fi — n)t and (Q — k) c are obtained from 
equation ( flUl) : 



3 ( r \ 3 

(^-^)t = ^K(r t )(^) 3 . (16) 

The growth rate of the wave- wave resonant instability depends on magnitudes of coupling terms 
and disk deformation (for a general expression for growth rate, see Kato et al. 2011). 
If the difference between Q and k is neglected, equation ( fl5i) gives 

"• = ^2^- (17) 



3 A rigorous way to know the relation among tl , r t and other quantities is to calculate the coupling terms and 
the resulting growth rate of oscillations for various sets of tl, r t and r c and to combine this procedure with 
that determining r t and r c in the following sections. This is, however, beyond the purpose of this paper. 
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Fig. 2. Schematic diagram showing frequencies and propagation regions of the u>\- and o^-oscillations. 
The scales of coordinates are arbitrary, and are not linear. The wi-oscillation is trapped between r c and 
r t . The inside of r c is the evanescent region. The a>2-oscillation has a negative frequency and propagates 
between the inner edge of the disk and tl- (Here, r\_, and rt are taken to be equal.) The propagation 
region is rather wide, but this does not mean that the a>2-oscillation is observed with large amplitude, 
since the oscillation has very short radial wavelength except near to r t . 

This is the results known as the parametric resonance, i.e., the 3 : 1 resonance in the case of 
?Tirj = 3, and the 2 : 1 resonance in the case of mrj = 4. It is noted that in the above cases of 
rriD = 3 and mo = 4, both of U2 and m>2 are negative. 

If we have a relation between r t and r c , equation ffl5l) gives r t and r c (and thus U\ 
and 0J2) as functions of q. The relation between r t and r c is obtained by considering that the 
wi-oscillation is trapped between r c and r t . That is, the relation is obtained by imposing the 
trapping condition, which is a relation among r c , r t and acoustic speed, c s . Thus, parameters 
to solve equation ( fl~5l) are q and c s . 

5. Trapping of One- Armed Low- Frequency Global Oscillation 

The one-armed low-frequency p-mode oscillation with frequency u\ is trapped in the 
region of r c < r < r t . The purpose here is to derive a relation between r c and r t by solving 
the eigenvalue problem, introducing acoustic speed in the disk. This problem has already been 
examined by Hirose & Osaki (1993) by using vertically integrated hydrodynamical equations. 
In different contexts, similar problems have been examined, i.e., by Okazaki (1991) for V/R 
variations of Be stars, and by Silbergleit et al. (1990) and others for c-mode oscillations 
in relativistic disks. Here, we consider the simplest situations that the radial wavelength of 
perturbations (on the order of r t — r c ) is so short that the characteristic radial scale of variations 
of unperturbed quantities in disks can be neglected except when we consider the radial variation 
of oj — (Q — k) . In this case the wave equation describing the trapped oscillation is (see Appendix) 
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o. 



where 



Q(r) 



(18) 



(19) 



Equation f lT8|) is now solved by the WKBJ method. In the region of r c < r < r t , the 
approximate solution can be written as 



Q- 1/4 exp ± i / Q 1/2 dr 



(20) 



The radius r c is the turning point of Q, where Q = 0. Near r c , the solution of equation (TPS!) 
thus can be expressed in terms of the Bessel functions of the order of ±1/3 (Morse & Feshbach 
1953). The asymptotic dependence of the solution for r ^> r c can be arranged so that it can 
be expressed in the form of equation (120~I) . Furthermore, by taking only the solution whose 
amplitude spatially damps in the region of r < r c (the evanescent region of oscillations), we 
have (Morse & Feshbach 1953) 



-1/4 



COS w 



12 



-7T 



COS 10 



12 



-7T 



where 



w(r)= / Q 1/2 dr. 



(21) 



(22) 



Next, we impose a boundary condition at r = r t . Since in this paper r t is taken to be 
the outer edge of the disks, the vanishing of the Lagrangian variation of pressure, i.e., 5p = 0, 
will be relevant. This is approximately equal to hi — (see Appendix for definition of hi), and 
thus to du r /dr = at r = r t . From the differentiation of the terms in the brackets of equation 
(f2~IT) with repect to r, we have (see also Silbergleit et al. 2001) 



w t 



n + -)■ 



(23) 



as the condition determining the wave trapping] 4 ! where n = 0, 1,2,3.... Since we are considering 
the fundamental mode in the radial direction, we adopt hereafter n — 0, and w t = 7r/4. 

We now perform the integration of Q 1 ^ 2 by using expression f|T9|) for Q. Considering 
that {oj — Q) 2 — k 2 ~ — 2Q[u — (Q — k)], and further that u = to\ = (Q — k) c and Q — k is given 
by equation (TTDjl . we have 



3(GM\ q 
V 2\DcVD* 



1- 



r c /D\ 3 / 2 
r/D ) 



(24) 



If we adopt u r = at r = r t as the boundary condition, we have 



w t 



Q^dr : 



n+ - r 
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Fig. 3. The q - r t /D and q - r c /D relations calculated by our model. The former relation is drawn for 
c s /(GM/Dy/ 2 = 0.02, which is almost free from the value of c s / (GM / D) 1 / 2 . The latter relation is drawn 
for three cases of c s /{GM/Df/ 2 = 0.01, 0.02, and 0.04. The region between r t and r c is the trapped region 
of the wi-oscillation. The trapped region becomes wide as c s increases. For comparison, the maximum 
radial size of periodic orbit of a test particle in the binary potential calculated by Paczhski (1977) is ploted 
for some values of q by black circles. The tidal instability occurs for q < 0.22. 

Here, we consider that (r t — r c )/r c <C 1 (the WKBJ method is still valid since GM/Dc 2 is a 
large quantity). Then, we can approximately perform the integration of equation (|22|) to have 



(GM \ V 2 (r t /Df/ 2 - (r c /Df' 2 



-<1 



which is changed to 

r t r c 2 w t 



(rjD) 1 / 2 



(25) 



(26) 



D D 3 q 1 / 2 (GM/D) 1 ! 2 ' 
This is a relation between r t /D and r c / D with dimentionless parameters c s / (GM/ D) l l 2 and q. 

6. Numerical Calculations 

Substitution of equations (fT6|) and (|26|) into equation (fT5l) gives the disk truncation 
radius r t as a function of q and c s / {GM / D) 1 / 2 . The q - r t relation for c s / {GM / D) 1 / 2 = 0.02 is 
shown in figure 3. The q - r t relation depends little on c s / (GM / D) 1 / 2 . The q - r c relation, on the 
other hand, depends on c s / (GM/D) 1 / 2 , and the relation in three cases of c s / (GM/D) 1 / 2 = 0.01, 
0.02, and 0.04 are also shown in figure 3. The region between r t and r c is the trapped one of 
the Wi-oscillation. As is expected, the width increases with increase of c s / '(GM / ' D) 1 ! 2 . 

Paczyhski (1977) estimated the maximum size of accretion disks in a binary system by 
calculating the maximum size of periodic orbit of a test particle. One of such radii, e.g., r max 
in his paper, is also shown in figure 3, for comparison, for some values of q. The comparison 
shows that the disk size is limitted by the tidal instability in the case where q < 0.22. 

The observed superhump frequency w s h is related to f2 or b and uj\ by (Osaki 1985) 
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Fig. 4. The q - e relation calculated for three cases of CsKGM/Df/ 2 = 0.01, 0.02, and 0.04. These curves 
are superposed on the diagram plotting observed q - e relation (hgure 34 of T. Kato et al. 2009). The 
names of nine sources in this figure are listed in table 6 of T.Kato et al. (2009). 



^sh — ^orb — (27) 

where f2 or b is the binary orbital frequency in the inertial frame [i.e., f2 or b = fi* rb (l + g) 1 ^ 2 ]. Then, 
the dimensionless quantity defined by e = P s h/-fobr _ 1, where P s h and P or b are the period of 
superhump and the orbital period (in the inertial frame), respectively, is written as 
_ c<;i/f2 or b 

1 — UJi/fl orh 

Here, in our present model, cui/f2 or b is given by [see eq. (ITT]) ] 
ui\ 3 q /r n \ 3 / 2 



(28) 



(29) 



n OTh A(i + q y/ 2 \Dj 

The q - e relation calculated by using the above relations (I28p and (T29]) is drawn in figure 4 for 
three cases of cJ^GM/D) 1 / 2 = 0.01, 0.02 and 0.04. In order to compare with observations, the 
curves are superposed on the q - e plot by T.Kato et al.(2012) for recently observed sources. 
The calculated curve in the case of c s / (GM / D) 1 / 2 = 0.02, which is c s /(_Df2* rb ), seems to be 
close to observational results. The dimensionless value of c s /(GMj D) l l 2 = 0.02 will be relevant 
observationally, since in accretion disks in cataclysmic variables, .Dfi orb ~ 500km s" 1 and c s ~ 
10km s" 1 . 

7. Discussions 

The present wave-wave resonant excitation model well describes the tidal instability - 
superhump phenomena in dwarf novae. This is not surprising, since the present excitation model 
can be regarded just as a combination of the tidal instability model by Lubow (1991a) and the 
precession of low frequency eccentric deformation by Osaki (1985) into a unified scheme. In the 
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tidal instability model by Lubow (1991a), the precession of the disk deformation is not essential 
in understanding the essence of the instability. The tidal torque and pressure force, however, 
make the deformation precess (Lubow 1992, 1994). Hence, there is a connection between 
the instability and eccentric precession, and much attention has been taken to the connection 
(Goodchild & Ogilvie 2006, Ogilvie 2007, Lubow 2010). In these studies, the temperature 
dependence of the precession has been examined (Lubow 2010). Here, we compare our results 
with those by Lubow (2010). Figure 5 shows the u\ - (H/r) t relation in the case of q = 0.1, 
where H is the disk half-thickness and taken simply to be c s /Qk and the subscript t denotes 
the value at the disk truncation radius, r t . The results show that the precession rate decreases 
with increase of (H/r) t , which agrees with results by Lubow (see figure 6 of Lubow 2010). 
Quantitatively, however, there are some differences. These differences partially come from 
differences of situations considered. In Lubow's calculations the disk edge is taken at 0.5D, 
while in our calculations it is taken at r t which is determined by relation ( TloT) and smaller 
than 0.5-D, depending on (H/r) t . In Lubow's calculations the rate of precession is taken to 
become the gravitational precession rate of a free particle on an excentric orbit in the limit of 
no temperature, while in our calculations it is taken to tend to (Q — n) t - 

It is noted that in our wave-wave resonant instability model two small amplitude oscil- 
lations are superposed on a steady deformed disk, and their resonant instability is considered. 
In other words, in a stage where the disk deformation is growing rapidly with time, the model 
cannot be applied. One of other important approximations involved in the present model is 
that r t and tl are assumed to be equal. In the limit of no temperature, the o^-oscillation is 
localized sharply at (the radius of the inner Lindblad resonance), and thus the coupling 
constants between the uj\- and 002- oscillations are large only when r t and tl are equal. In 
disks with finite temperature, however, the radial region where the coupling constants are large 
becomes wide. As is mentioned in a footnote, to know a proper relation between r t and 7l we 
must calculate the radial distribution of the coupling coefficients, using the functional forms 
of the oscillations. This is a complicated problem. This situation corresponds to the fact that 
in disks with finite temperature the resonant region has a finite width (e.g., Meyer- Vernet & 
Sicardy 1987). 

In our model the presence of two oscillation modes with opposite signs of oj — mfl are 
essential for instability, and one of the oscillations excited (i.e., one-armed low-frequency oscilla- 
tion with bj — mQ < 0) brings about a precession of the deformation. The other oscillation with 
uj — mQ > is not directly related to the superhump, but important for resonant instability. 
This latter mode corresponds to the intermediate disk deformation mode in Lubow's resonant 
feedback process (Lubow 1991a). In our present model, |m 2 | = 2 like Lubow's case, and its 
frequency j^l changes with the change of u\, since they are related by a resonance condition 
as \oj 2 \ =u>£) — u)i. The frequency \u 2 \ of the oscillation is shown in figure 6 as a function of q 
in the case of c s / (GM/D) 1 / 2 = 0.02. The curve depends little on the value of c s / \GM '/ ' D) 1 ' 2 , 
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Fig. 5. The temperature dependence of u>\ in the case of q = 0.1. The rate of precession, uj\, is normalized 
by fi or b and the effects of temperature are measured by (H/r)t, where H is the disk half-thickness and 
H/r is taken to be c s /£!k- In the limit of zero temperature, u>\ tends to (f2 — n) t in our model. For 
comparison, fi-K at the radius of = 3f2* bs is shown by a thin line. 

but changes in the range of 3.0 ~ 2.4f2 or b for change of q in the range of < q < 0.4. 

Observational detection of the w 2 -oscillation and its q-dependence is an interesting sub- 
ject to judge whether the present wave-wave resonant instability model represents real situations 
of tidal instability - superhump phenomena. Osaki (2003) pointed out that complex superhump 
light curves for the 2001 outburst of WZ Sagittae can be interpreted as excitation of many oscil- 
lation modes. Especially, he emphasizes the appearance of two armed oscillation with frequency 
3f2 or b — uji, and considers it as a support of Lubow' mode-coupling model. Except for this case, 
however, there seems to be no observational evidence which suggests the coexistence (with the 
Wi-oscillation) of the oj 2 -oscillation in the superoutburst stage (private communication by T. 
Kato). Osaki suggests that this is due to the fact that many sources with superhumps are 
observed pole-on. In addition there will be following situations for difficulty of observations 
of the (^-oscillation. As mentioned before, the ^-oscillation has short wavelength except in 
the region near to r t , although its propagation region is wide. Hence, a large luminosity vari- 
ation will not be expected, since various phases of the oscillation are superposed on observed 
quantities. Furthermore, r*L = r t adopted in this paper is an approximation. The outer part of 
disks will be fluctuating with time in general, and this will introduce time fluctuations in the 
relation between tl and r t . If so, the phase and frequency of the ^-oscillation will be subject 
to it, more than those of the o;i-oscillation. This will smear observational properties of the 
^-oscillation. 

So far, we have considered the so-called 3 : 1 resonance which occurs for = 3. In 
this case the cu2-oscillation is a two-armed one (i.e., |rri2 1 = 2) with frequency close to 3f2 or b as 
mentioned above. For tidal deformation of mo = 4, on the other hand, the instability occurs 
near 2 : 1 resonance [see equation (fTTj) ]. In this case, r t is larger than that in the case of the 
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Fig. 6. The q - \lj 2 \ relation calculated for c s j (GM j D) 1 ! 2 = 0.02. The corresponding curves in the cases 
of c s /(GM /D) 1 / 2 = 0.01, and 0.04 are almost the same as that in the case of c s / (GM/D) 1 / 2 = 0.02. 
For reference, the frequency |oj2 1 i n the case of the 2 : 1 resonance is also drawn as a function of q for 
c s /{GM/D) 1 f 2 = Qm. 

3 : 1 resonance. The ^-oscillation in this case is three-armed (i.e., \m,2\ = 3) and its frequency 
is close to 4f2 orb [see equation ffT4"l) ]. As mentioned before, observational detection of the u 2 - 
oscillation is interesting to evaluate validity of the present wave-wave resonant model. Thus, 
for reference, the frequency, \u 2 \, in the case of the 2 : 1 resonance is also drawn in figure 5 as 
a function of q. 

We think that the high frequency quasi-periodic oscillations (HF QPOs) in microquasars 
will result from a similar process considered here. Hence, detailed examinations how well the 
present wave-wave resonant model can describe the superoutburst-superhump phenomena in 
dwarf novae is helpful to consider more the excitation mechanism of HF QPOs in BH LMXBs. 

The author thanks T. Kato for discussions on various observational aspects of super- 
humps, and Y. Osaki for pointing out the presence of high frequency oscillation of f2 orb — u>i in 
an outburst stage of WZ Sagittae. The author also thanks the referee for pointing out recent 
theoretical developments on disk precession and for many helpful suggestions. 

Appendix: Wave Equation Describing Disk Oscillations 

For simplicity, the disk is assumed to be vertically isothermal. Then, the vertical hydro- 
static balance gives 



where p (r,z) is the density in the unperturbed disk and Poo( r ) is that on the equatorial plane. 
The scale height in the vertical direction, H, is given by 




(30) 




(31) 
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i(u - mVt)pi + -r-(rpoMr) - i—pou v 



where c s (r) is the isothermal acoustic speed and Qj_ is the vertical epicyclic frequency in the 
tidally deformed disks. 

On the disks we impose small-amplitude isothermal perturbations. The perturbations 
are assumed to be proportional to exp[i(ut — rrup)]. Then, hydrodynamical equations describing 
the small-amplitude perturbations are 

d , m 

IT 1 
ror 

+J^(A>u*) = 0, (32) 
dh\ 

i(oj — mfl)u r — 2Qu v = — — , (33) 

Or 

2 

/ . /tj TTl . > 

i(u) — m\l)u v + — u r = i — hi, (34) 
i{u — mVt)u z = — -77— , (35) 



dz 



where h\ is defined by 



h 1 = ^ = c 2 A. (36) 
Po Po 

Here, ( velocity perturbations over unperturbed one (0, rfi, 0), and pi and pi 

are the density and pressure perturbations over the unperturbed ones, po and po, respectively. 

Hereafter, we consider perturbations whose radial wavelength is shorter than the char- 
acteristic radial scales of unperturbed quantities in disks. Elimination of from equation (1331) 
and (1341) then gives 

dh 

{-(u -mtt) 2 + K 2 ]u r = -i(cu -mtt)— -. (37) 

or 



Furthermore, the continuity equation ( 1321) is reduced to 

hi du r ( d z \ , , 

''"-"•^^ns-ff)"- - (38) 

Under the same approximation, elimination of u z from this equation by using equation (135*1) 
gives 

d 2 hi z dhi (u-mQ) 2 du r 

^-Jp^ + ~^— hl = t(U ~ mn) ^- (39) 
Equations ( 1371) and ( 1391) are the set of equations to be solved. As mentioned in the text, 
we are interested in one-armed (m — 1) low- frequency oscillations, where u ~ Q — k. Hence, 
the radial variation of (ou — Q) 2 — k 2 should be taken into account, but the radial variations 
of other quantities are neglected. Then, operating d 2 /dz 2 — (z/H 2 )d/dz + (u — mQ) 2 /c 2 on 
equation (|37j) and eliminating hi from the resulting equation by using (139|) . we have an partial 
differential equation with respect to u r . After manipulations we have 

d 2 z d\ (u-ny 
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dz 2 H 2 dz J Ur c? Ur 



;u r = 0, (40) 



(u — f2) 2 — k 2 dr 2 
where m = 1 has been adopted. 

This partial differential equation can be easily decomposed into two ordinary differential 
equations. That is, decomposing u r (r,z) as u r (r,z) = g(z)f(r), we have 

^-4£ + *W) = (41) 



, dz 2 H 2 dz 
and 

(oo-n) 2 d 2 (u-Q) 2 



(uj — Vt) 2 — k 2 dr 2 c. 



K 



fir) = 0, (42) 



where K is the separation constant, and can be determined by solving equation ( l4Tj) with 
relevant boundary conditions at z = ±oo. 

Equation ( I4T1) is the Hermite equation. Imposing that p (z)g (which is proportional 
to kinetic energy of perturbations) remains finite at z = ±oo, we find that K should be zero or 
a positive integer and g(z) is the Hermite polynomials (Okazaki et al. 1987): 

K = n, n = 0,l,2,... (43) 

and 

g(z)=H n (z/H). (44) 

In the present problem, we are interested in the fundamental p-mode oscillation, i.e., 
n = 0. Hence, equation fT42|) is reduced to 
T d 2 (cu-n) 2 -K 2 ] , . 

^~^h — ] /(r) = a (45) 

This is the equation which we treat in the text. 
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Note added on Jan. 3, 2013 

In the present paper, we have restricted our attention only to the resonances of w 2 = (jJi ± 
and m 2 = m 1 ± m^. If the resonances of u = —U\ ± u-^ and m 2 = —vn\ ± are considered, 
the instability condition in the latter cases is found to be (-&l/o;i)(-E 2 /c<; 2 ) > 0, different from 
(Ei/ui)(E 2 /uj2) < in the former cases of w 2 = u)\ ± c^d and m 2 = mi ±m^. If the latter 
resonances are taken into account, we see that in addition to the (^-oscillations considered in 
this paper (i.e., u 2 = u± — wd < and m 2 = mi — mo < 0), the w 2 -mode of oscillations with 
uj 2 = —tui + ud > and m 2 = —mi + mo > is also excited. These two oscillations are, however, 
the same, since both signs of u and m are changed simultaneously. A detailed re-examination 
of wave-wave resonant instability by Kato et al. (2011) will be made in a subsequent paper. 
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